Expression profile of the entire detoxification gene inventory of the western honeybee, Apis mellifera across life stages.

The western honeybee, Apis mellifera, is a managed pollinator of many crops and potentially exposed to a wide range of foreign compounds, including pesticides throughout its life cycle. Honeybees as well as other insects recruit molecular defense mechanisms to facilitate the detoxification of xenobiotic compounds. The inventory of detoxification genes (DETOXome) is comprised of five protein superfamilies: cytochrome P450 monooxygenases (P450), carboxylesterases, glutathione S-transferases (GST), UDP-glycosyl transferases (UGT) and ATP-binding cassette (ABC) transporters. Here we characterized the gene expression profile of the entire honeybee DETOXome by analyzing 47 transcriptomes across the honeybee life cycle, including different larval instars, pupae, and adults. All life stages were well separated by principal component analysis, and K-means clustering revealed distinct temporal patterns of gene expression. Indeed, >50% of the honeybee detoxification gene inventory is found in one cluster and follows strikingly similar expression profiles, i.e., increased expression during larval development, followed by a sharp decline after pupation and a steep increase again in adults. This cluster includes 29 P450 genes dominated by CYP3 and CYP4 clan members, 15 ABC transporter genes mostly belonging to the ABCC subfamily and 13 carboxylesterase genes including almost all members involved in dietary/detox and hormone/semiochemical processing. RT-qPCR analysis of selected detoxification genes from all families revealed high expression levels in various tissues, especially Malpighian tubules, fatbody and midgut, supporting the view that these tissues are essential for metabolic clearance of environmental toxins and pollutants in honeybees. Our study is meant to spark further research on the molecular basis of detoxification in this critical pollinator to better understand and evaluate negative impacts from potentially toxic substances. Additionally, the entire gene set of 47 transcriptomes collected and analyzed provides a valuable resource for future honeybee research across different disciplines.

A conserved hymenopteran-specific family of cytochrome P450s protects bee pollinators from toxic nectar alkaloids.

Many plants produce chemical defense compounds as protection against antagonistic herbivores. However, how beneficial insects such as pollinators deal with the presence of these potentially toxic chemicals in nectar and pollen is poorly understood. Here, we characterize a conserved mechanism of plant secondary metabolite detoxification in the Hymenoptera, an order that contains numerous highly beneficial insects. Using phylogenetic and functional approaches, we show that the CYP336 family of cytochrome P450 enzymes detoxifies alkaloids, a group of potent natural insecticides, in honeybees and other hymenopteran species that diverged over 281 million years. We linked this function to an aspartic acid residue within the main access channel of CYP336 enzymes that is highly conserved within this P450 family. Together, these results provide detailed insights into the evolution of P450s as a key component of detoxification systems in hymenopteran species and reveal the molecular basis of adaptations arising from interactions between plants and beneficial insects.

The mustard leaf beetle, Phaedon cochleariae, as a screening model for exogenous RNAi-based control of coleopteran pests.

RNA interference (RNAi) is a promising, selective pest control technology based on the silencing of targeted genes mediated by the degradation of mRNA after the ingestion of double-stranded (ds) RNA. However, the identification of the best target genes remains a challenge, because large scale screening is only feasible in lab model systems and it remains unclear, to what degree such data can be transferred to pest species. Here, we report on our efforts to transfer target genes found in a lab model to the mustard leaf beetle, Phaedon cochleariae. The mustard leaf beetle can be reared easily and resource-efficient in large quantities all year round and is an established chrysomelid pest for higher throughput screening approaches in the crop protection industry. Mustard leaf beetle transcriptome sequencing and assembly revealed genes orthologous to those previously described as highly efficient RNAi targets in the model beetle Tribolium castaneum. First, we observed mortality after injection of dsRNA targeting the respective orthologous genes in 2nd instar mustard beetle larvae. Next, we adopted a robust, automated multi-well plate foliar RNAi screening procedure with 2nd instar larvae of the mustard leaf beetle to assess those genes. Indeed, foliar application and oral uptake of dsRNA targeting the same genes resulted in larval mortality as well. The most effective target genes with a strong (lethal) phenotype - at dsRNA doses as low as 300 ng/leaf disc (equal to 9.6 g/ha) - were srp54k, rop, αSNAP, rpn7 and rpt3. Rather limited effects were observed after application of dsRNA targeting cactus, shibire and PP-α, though they had previously been shown to be highly lethal in red flour beetle. Importantly, our experiments demonstrated that the overall efficacy pattern obtained after oral dsRNA application was well correlated with the results obtained after dsRNA injection. RT-qPCR confirmed significant target gene knock-down after normalization by employing three reference genes shown to be stably expressed across life stages. In summary, several RNAi targeted genes elicited a strong lethal phenotype and significant target gene knock-down after feeding, suggesting P. cochleariae as a potential coleopteran screening model for foliarly applied exogenous RNAi.

A toxicogenomics approach reveals characteristics supporting the honey bee (Apis mellifera L.) safety profile of the butenolide insecticide flupyradifurone.

Flupyradifurone, a novel butenolide insecticide, selectively targets insect nicotinic acetylcholine receptors (nAChRs), comparable to structurally different insecticidal chemotypes such as neonicotinoids and sulfoximines. However, flupyradifurone was shown in acute toxicity tests to be several orders of magnitude less toxic to western honey bee (Apis mellifera L.) than many other insecticides targeting insect nAChRs. The underlying reasons for this difference in toxicity remains unknown and were investigated here. Pharmacokinetic studies after contact application of [14C]flupyradifurone to honey bees revealed slow uptake, with internalized compound degraded into a few metabolites that are all practically non-toxic to honey bees in both oral and contact bioassays. Furthermore, receptor binding studies revealed a lack of high-affinity binding of these metabolites to honey bee nAChRs. Screening of a library of 27 heterologously expressed honey bee cytochrome P450 enzymes (P450s) identified three P450s involved in the detoxification of flupyradifurone: CYP6AQ1, CYP9Q2 and CYP9Q3. Transgenic Drosophila lines ectopically expressing CYP9Q2 and CYP9Q3 were significantly less susceptible to flupyradifurone when compared to control flies, confirming the importance of these P450s for flupyradifurone metabolism in honey bees. Biochemical assays using the fluorescent probe substrate 7-benzyloxymethoxy-4-(trifluoromethyl)-coumarin (BOMFC) indicated a weak, non-competitive inhibition of BOMFC metabolism by flupyradifurone. In contrast, the azole fungicides prochloraz and propiconazole were strong nanomolar inhibitors of these flupyradifurone metabolizing P450s, explaining their highly synergistic effects in combination with flupyradifurone as demonstrated in acute laboratory contact toxicity tests of adult bees. Interestingly, the azole fungicide prothioconazole is only slightly synergistic in combination with flupyradifurone - an observation supported by molecular P450 inhibition assays. Such molecular assays have value in the prediction of potential risks posed to bees by flupyradifurone mixture partners under applied conditions. Quantitative PCR confirmed the expression of the identified P450 genes in all honey bee life-stages, with highest expression levels observed in late larvae and adults, suggesting honey bees have the capacity to metabolize flupyradifurone across all life-stages. These findings provide a biochemical explanation for the low intrinsic toxicity of flupyradifurone to honey bees and offer a new, more holistic approach to support bee pollinator risk assessment by molecular means.

The genetic architecture of a host shift: An adaptive walk protected an aphid and its endosymbiont from plant chemical defenses.

Host shifts can lead to ecological speciation and the emergence of new pests and pathogens. However, the mutational events that facilitate the exploitation of novel hosts are poorly understood. Here, we characterize an adaptive walk underpinning the host shift of the aphid Myzus persicae to tobacco, including evolution of mechanisms that overcame tobacco chemical defenses. A series of mutational events added as many as 1.5 million nucleotides to the genome of the tobacco-adapted subspecies, M. p. nicotianae, and yielded profound increases in expression of an enzyme that efficiently detoxifies nicotine, both in aphid gut tissue and in the bacteriocytes housing the obligate aphid symbiont Buchnera aphidicola. This dual evolutionary solution overcame the challenge of preserving fitness of a mutualistic symbiosis during adaptation to a toxic novel host. Our results reveal the intricate processes by which genetic novelty can arise and drive the evolution of key innovations required for ecological adaptation.

Identification and functional characterization of a novel acetyl-CoA carboxylase mutation associated with ketoenol resistance in Bemisia tabaci.

Insecticides of the tetronic/tetramic acid family (cyclic ketoenols) are widely used to control sucking pests such as whiteflies, aphids and mites. They act as inhibitors of acetyl-CoA carboxylase (ACC), a key enzyme for lipid biosynthesis across taxa. While it is well documented that plant ACCs targeted by herbicides have developed resistance associated with mutations at the carboxyltransferase (CT) domain, resistance to ketoenols in invertebrate pests has been previously associated either with metabolic resistance or with non-validated candidate mutations in different ACC domains. A recent study revealed high levels of spiromesifen and spirotetramat resistance in Spanish field populations of the whitefly Bemisia tabaci that was not thought to be associated with metabolic resistance. We confirm the presence of high resistance levels (up to >640-fold) against ketoenol insecticides in both Spanish and Australian B. tabaci strains of the MED and MEAM1 species, respectively. RNAseq analysis revealed the presence of an ACC variant bearing a mutation that results in an amino acid substitution, A2083V, in a highly conserved region of the CT domain. F1 progeny resulting from reciprocal crosses between susceptible and resistant lines are almost fully resistant, suggesting an autosomal dominant mode of inheritance. In order to functionally investigate the contribution of this mutation and other candidate mutations previously reported in resistance phenotypes, we used CRISPR/Cas9 to generate genome modified Drosophila lines. Toxicity bioassays using multiple transgenic fly lines confirmed that A2083V causes high levels of resistance to commercial ketoenols. We therefore developed a pyrosequencing-based diagnostic assay to map the spread of the resistance alleles in field-collected samples from Spain. Our screening confirmed the presence of target-site resistance in numerous field-populations collected in Sevilla, Murcia and Almeria. This emphasizes the importance of implementing appropriate resistance management strategies to prevent or slow the spread of resistance through global whitefly populations.

Molecular characterization of Cry1F resistance in fall armyworm, Spodoptera frugiperda from Brazil.

Fall armyworm, Spodoptera frugiperda (J.E. Smith) is a major lepidopteran pest of maize in Brazil and its control particularly relies on the use of genetically engineered crops expressing Bacillus thuringiensis (Bt) toxins such as Cry1F. However, control failures compromising the efficacy of this technology have been reported in many regions in Brazil, but the mechanism of Cry1F resistance in Brazilian fall armyworm populations remained elusive. Here we investigated the molecular mechanism of Cry1F resistance in two field-collected strains of S. frugiperda from Brazil exhibiting high levels of Cry1F resistance. We first rigorously evaluated several candidate reference genes for normalization of gene expression data across strains, larval instars and gut tissues, and identified ribosomal proteins L10, L17 and RPS3A to be most suitable. We then investigated the expression pattern of ten potential Bt toxin receptors/enzymes in both neonates and 2nd instar gut tissue of Cry1F resistant fall armyworm strains compared to a susceptible strain. Next we sequenced the ATP-dependent Binding Cassette subfamily C2 gene (ABCC2) and identified three mutated sites present in ABCC2 of both Cry1F resistant strains: two of them, a GY deletion (positions 788-789) and a P799 K/R amino acid substitution, located in a conserved region of ABCC2 extracellular loop 4 (EC4) and another amino acid substitution, G1088D, but in a less conserved region. We further characterized the role of the novel mutations present in EC4 by functionally expressing both wild type and mutated ABCC2 transporters in insect cell lines, and confirmed a critical role of both sites for Cry1F binding by cell viability assays. Finally, we assessed the frequency of the mutant alleles by pooled population sequencing and pyrosequencing in 40 fall armyworm populations collected from maize fields in different regions in Brazil. We found that the GY deletion being present at high frequency. However we also observed many rare alleles which disrupt residues between sites 783-799, and their diversity and abundance in field collected populations lends further support to the importance of the EC4 domain for Cry1F toxicity.

Unravelling the Molecular Determinants of Bee Sensitivity to Neonicotinoid Insecticides.

The impact of neonicotinoid insecticides on the health of bee pollinators is a topic of intensive research and considerable current debate [1]. As insecticides, certain neonicotinoids, i.e., N-nitroguanidine compounds such as imidacloprid and thiamethoxam, are as intrinsically toxic to bees as to the insect pests they target. However, this is not the case for all neonicotinoids, with honeybees orders of magnitude less sensitive to N-cyanoamidine compounds such as thiacloprid [2]. Although previous work has suggested that this is due to rapid metabolism of these compounds [2-5], the specific gene(s) or enzyme(s) involved remain unknown. Here, we show that the sensitivity of the two most economically important bee species to neonicotinoids is determined by cytochrome P450s of the CYP9Q subfamily. Radioligand binding and inhibitor assays showed that variation in honeybee sensitivity to N-nitroguanidine and N-cyanoamidine neonicotinoids does not reside in differences in their affinity for the receptor but rather in divergent metabolism by P450s. Functional expression of the entire CYP3 clade of P450s from honeybees identified a single P450, CYP9Q3, that metabolizes thiacloprid with high efficiency but has little activity against imidacloprid. We demonstrate that bumble bees also exhibit profound differences in their sensitivity to different neonicotinoids, and we identify CYP9Q4 as a functional ortholog of honeybee CYP9Q3 and a key metabolic determinant of neonicotinoid sensitivity in this species. Our results demonstrate that bee pollinators are equipped with biochemical defense systems that define their sensitivity to insecticides and this knowledge can be leveraged to safeguard bee health.

Neofunctionalization of Duplicated P450 Genes Drives the Evolution of Insecticide Resistance in the Brown Planthopper.

Gene duplication is a major source of genetic variation that has been shown to underpin the evolution of a wide range of adaptive traits [1, 2]. For example, duplication or amplification of genes encoding detoxification enzymes has been shown to play an important role in the evolution of insecticide resistance [3-5]. In this context, gene duplication performs an adaptive function as a result of its effects on gene dosage and not as a source of functional novelty [3, 6-8]. Here, we show that duplication and neofunctionalization of a cytochrome P450, CYP6ER1, led to the evolution of insecticide resistance in the brown planthopper. Considerable genetic variation was observed in the coding sequence of CYP6ER1 in populations of brown planthopper collected from across Asia, but just two sequence variants are highly overexpressed in resistant strains and metabolize imidacloprid. Both variants are characterized by profound amino-acid alterations in substrate recognition sites, and the introduction of these mutations into a susceptible P450 sequence is sufficient to confer resistance. CYP6ER1 is duplicated in resistant strains with individuals carrying paralogs with and without the gain-of-function mutations. Despite numerical parity in the genome, the susceptible and mutant copies exhibit marked asymmetry in their expression with the resistant paralogs overexpressed. In the primary resistance-conferring CYP6ER1 variant, this results from an extended region of novel sequence upstream of the gene that provides enhanced expression. Our findings illustrate the versatility of gene duplication in providing opportunities for functional and regulatory innovation during the evolution of an adaptive trait.

Development of a lateral flow test to detect metabolic resistance in Bemisia tabaci mediated by CYP6CM1, a cytochrome P450 with broad spectrum catalytic efficiency.

Cotton whitefly, Bemisia tabaci (Genn.) (Homoptera: Aleyrodidae) is a major sucking pest in many agricultural and horticultural cropping systems globally. The frequent use of insecticides of different mode of action classes resulted in populations resisting treatments used to keep numbers under economic damage thresholds. Recently it was shown that resistance to neonicotinoids such as imidacloprid is linked to the over-expression of CYP6CM1, a cytochrome P450 monooxygenase detoxifying imidacloprid and other neonicotinoid insecticides when recombinantly expressed in insect cells. However over-expression of CYP6CM1 is also known to confer cross-resistance to pymetrozine, an insecticide not belonging to the chemical class of neonicotinoids. In addition we were able to demonstrate by LC-MS/MS analysis the metabolisation of pyriproxyfen by recombinantly expressed CYP6CM1. Based on our results CYP6CM1 is one of the most versatile detoxification enzymes yet identified in a pest of agricultural importance, as it detoxifies a diverse range of chemical classes used to control whiteflies. Therefore we developed a field-diagnostic antibody-based lateral flow assay which detects CYP6CM1 protein at levels providing resistance to neonicotinoids and other insecticides. The ELISA based test kit can be used as a diagnostic tool to support resistance management strategies based on the alternation of different modes of action of insecticides.

Over-expression of cytochrome P450 CYP6CM1 is associated with high resistance to imidacloprid in the B and Q biotypes of Bemisia tabaci (Hemiptera: Aleyrodidae).

The two most damaging biotypes of Bemisia tabaci, B and Q, have both evolved strong resistance to the neonicotinoid insecticide imidacloprid. The major mechanism in all samples investigated so far appeared to be enhanced detoxification by cytochrome P450s monooxygenases (P450s). In this study, a polymerase chain reaction (PCR) technology using degenerate primers based on conserved P450 helix I and heme-binding regions was employed to identify P450 cDNA sequences in B. tabaci that might be involved in imidacloprid resistance. Eleven distinct P450 cDNA sequences were isolated and classified as members of the CYP4 or CYP6 families. The mRNA expression levels of all 11 genes were compared by real-time quantitative RT-PCR across nine B and Q field-derived strains of B. tabaci showing strong resistance, moderate resistance or susceptibility to imidacloprid. We found that constitutive over-expression (up to approximately 17-fold) of a single P450 gene, CYP6CM1, was tightly related to imidacloprid resistance in both the B and Q biotypes. Next, we identified three single-nucleotide polymorphic (SNP) markers in the intron region of CYP6CM1 that discriminate between the resistant and susceptible Q-biotype CYP6CM1 alleles (r-Q and s-Q, respectively), and used a heterogeneous strain to test for association between r-Q and resistance. While survivors of a low imidacloprid dose carried both the r-Q and s-Q alleles, approximately 95% of the survivors of a high imidacloprid dose carried only the r-Q allele. Together with previous evidence, the results reported here identify enhanced activity of P450s as the major mechanism of imidacloprid resistance in B. tabaci, and the CYP6CM1 gene as a leading target for DNA-based screening for resistance to imidacloprid and possibly other neonicotinoids in field populations.

Multiple origins of pyrethroid resistance in sympatric biotypes of Bemisia tabaci (Hemiptera: Aleyrodidae).

The two most damaging biotypes of Bemisia tabaci, B and Q, are sympatric in the Mediterranean basin and show high resistance to pyrethroids synergized by organophosphates. Previous work showed that in the B biotype, this resistance is associated with the L925I mutation in the para-type voltage gated sodium channel. Here we identified two mutations in the para-type voltage gated sodium channel associated with resistance to pyrethroids synergized by organophosphates in the Q biotype: the L925I mutation that occurs in the B biotype, and substitution of threonine to valine, at position 929 (T929V). To determine if the L925I and T929V mutations have single or multiple origins, we sequenced the DNA regions flanking the mutations from 13 B and Q strains collected worldwide. The survey identified five resistant alleles and five susceptible alleles. In the resistant alleles, the nucleotide diversity was low within biotypes (0.001), but high between biotypes (0.033). Nucleotide diversity in susceptible alleles was high between the two biotypes (0.028). These observations are consistent with multiple independent origins of resistance. Although the B and Q biotypes coexist in several regions of the Mediterranean basin, divergence in their DNA sequences at the para-type voltage gated sodium channel locus suggests gene flow between these biotypes is low or nil.